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Abstract. The Standard Model of particle physics and the theory of General 
Relativity (GR) currently provide a good description of almost all phenomena of 
particle physics and gravitation that have received controlled experimental tests. 
However, the Standard Model contains many a priori variable parameters whose 
values, and whose apparent (near-)constancy, have yet to receive a convincing 
theoretical explanation. At the same time, GR may now require to be extended or 
altered at the largest length scales, to account for the recent apparent accelerated 
cosmological expansion. In this introductory review I present theoretical aspects 
of the search for explanations of the values and possible variations of fundamental 
"constants", focusing on the possibility of unification of interactions. I also relate 
cosmological variations to modifications of gravity both locally and cosmologically. 

1 Introduction 

1.1 Motivation: Cosmological variation has occurred! 

It is a fundamental postulate of General Relativity (GR) that the same "non-gravitational" , 
local experiment will give the same result at all points in space-time. This is known as the prin- 
ciple of Local Position Invariance (LPI) , and forms part of the Einstein Equivalence Principle 
(EEP)[^]Yet it is very likely that during the previous evolution of the Universe, elementary 
particle properties, in particular masses and couplings, were substantially different from their 
present values. In a hot Big Bang at a temperature above about 1 TeV, it is believed that 
the Higgs field undergoes a phase transition, such that its expectation value becomes zero (as 
opposed to its present-day value of approximately 247 GeV) and thus the leading contribution 
to most elementary particle masses would vanish. Particle physics under such conditions would 
be utterly different from today. 

This is in one sense an artificial example, since the variation in the expectation value of 
the Higgs is triggered by temperature, whereas the LPI is supposed to hold after extraneous 
factors such as temperature are removed. Yet it does show that in the presence of scalar fields 
that may take different values at different points in spacetime, the assumptions underlying 
current theories of particle physics and gravitation may become invalid. It is also a warning 
that properties that are locally valid to high accuracy (i.e. constancy of particle properties) 
need not be so when extrapolated over very large distances or cosmological time intervals. 
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1.2 Fundamental parameters in particle physics and gravity 

The modern understanding of elementary particle physics and gravitation is encoded in the 
quantum field theoretical construction of the Standard Model, and in Einstein's theory of GR 
respectively. The Standard Model has one dimensionful parameter, the Higgs v.e.v. (<fi), and 
on the order of 20 dimensionless parameters which describe the static properties of particles 
and their interactions via gauge (vector) boson exchange. Gravity introduces one further di- 
mensionful parameter, the Newton constant, which may be exchanged for the Planck mass 
Mp = (SttGn^ 1 / 2 (in "natural" units where h and c are unity). Since the value of a dimen- 
sionful parameter is only known in terms of units, which themselves have to be defined by some 
physical system (for example a standard length or mass), the fundamental parameters whose 
values we can meaningfully discuss include only the ratio of the Higgs v.e.v. to the Planck mass. 

The Standard Model is believed to describe particle interactions from the smallest energies 
(corresponding to the largest length scales: only electromagnetism can act over macroscopic 
distances) up to energies at least of order 100 GeV, and possibly much higher. Whether new 
physics will be needed to supplement or supplant the SM at higher energies is an experimental 
question that is currently being investigated at the TeVatron at Fermilab and which the forth- 
coming operation of the LHC at CERN is intended to address directly. If new physics appears 
then the parameters of the SM may be expressed in terms of those of the new theory: this may 
lead to an increase, or (hopefully, from the point of view of explanatory power) a decrease in 
the number of independent parameters. 

One paradigm of how this may happen is unification: the variegation of different particles 
and interactions may be a symptom of the energy scale we are capable of observing them at. 
There may be a larger symmetry than we currently observe, which is realized in a nontrivial 
way or effectively hidden, but would be manifest if we were able to observe interactions at much 
higher energies. The consequences for fundamental parameters are explained in Section[3j under 
simple unified models we would expect variations to occur in all parameters in a correlated 
fashion. Hence a nonzero variation could test and distinguish between unified models if the 
observations constrain more than one fundamental parameter. 

1.3 Gravitational and cosmological related effects 

I return to the gravitational sector in Section |4j The direct physical motivation for consider- 
ing gravitational effects of spacetime variation is that if the forces that hold matter together 
depend on position, bodies will be accelerated towards regions where the mass-energy of their 
internal structure is lower, in addition to the GR gravitational acceleration. This is interesting, 
and potentially problematic, since these anomalous accelerations break the strong equivalence 
principle (which tells us that the only gravitational force is due to GR and that the strength 
of gravity is everywhere constant) and also in general the weak equivalence principle (which 
tells us that gravitational forces produce the same acceleration of a test body regardless of its 
composition) . 

To get a more quantitative view of equivalence-principle violations, physical models of space- 
time variation are required. Within a Lorentz-covariant framework this is only possible by 
introducing a new degree of freedom — usually a scalar field — that sources the variation. The 
choice of model determines the action of this degree of freedom and its couplings to gravity, 
the Standard Model and possibly other sectors (such as dark matter). Some of the strongest 
bounds on models of variation then arise from the effects of scalar-mediated forces. Indeed, if a 
scalar is to source a variation over time scales of the order of the age of the Universe, it should 
be extremely light and thus mediate very long-range interactions. 

Such new degrees of freedom also have nontrivial consequences for the global behaviour 
of cosmology, since they may contribute nontrivially to the potential and kinetic energy (i.e. 
stress-energy tensor) at cosmologically large scales. This may either be used as a candidate 
explanation for the apparent accelerated expansion at recent epochs, as so-called dark energy, 
or to bound the possible behaviour of the scalar over cosmic time, in particular to set an 
observational upper limit on its rate of change. 



Will be inserted by the editor 



3 



Given particular models of cosmological variation, one can test them by looking at epochs 
in cosmology where the variation is expected to be large, or where its effects may be accessible 
to observation. The earliest time at which bounds can be put on the values of particle physics 
parameters is that of primordial nucleosynthesis (BBN) , which depends on all three fundamental 
forces (strong, electromagnetic and weak) and on gravity, via the expansion rate of the Universe. 
Although such bounds are generally at the level of 10 -2 -10 -3 for fractional variations pQ thus 
not directly competitive with later spectroscopic results, they apply at a redshift z ~ 10 10 
where the variation might be considerably larger. 



2 The Standard Model and its parameters 

The definition of the Standard Model comprises the enumeration of matter particles taken to 
be pointlikc fcrmions, and the description of their interactions via bosons, with the appropriate 
coupling strengths. These interactions are divided into electromagnetic, weak and strong forces: 
the essential feature of the Standard Model is that the electromagnetic and weak interactions 
are both relics of a more fundamental set of forces, technically an SU(2)xU(l) gauge theory, 
whose effects are masked by the large nonzero value of the Higgs v.e.v., as will be described 
shortly. 

The strong force also results from a gauge theory, quantum chromodynamics (QCD): its ele- 
mentary excitations, gluons and quarks, are however not individually observable, being confined 
by the strong nonlinear self-interactions of the gluons. The properties of strongly interacting 
particles (hadrons) are then believed to be explained as "colour-neutral" bound states of quarks 
and gluons; nonperturbative techniques in quantum field theory such as the lattice are required 
to calculate them. Evidence for individual quarks and gluons only appears at high energies 



(well above the mass scales of bound states) where, as explained in Section 3.1 the effective 
coupling strength becomes weak. 

The elementary particle content of the SM can be written via the SU(3)c><SU(2)vi/xU(l)y 
quantum numbers: 

3 x {q(3, 2) 1/6 + [ u c (3, l)_ 2/3 + d c (3, l) 1/3 ] + 1(1, 2)_ 1/2 + [e c (l, 1)! + 1)„]} (1) 

where each representation is written as a left-handed Weyl fermion and we are including the 
possible presence of a "right-handed neutrino" i.e. a neutral, SU(2) singlet Weyl fermion. The 
factors of 3 simply denote the three generations, which appear to be identical apart from 
their mass terms. As required when electromagnetism is a mixture of the SU(2) and U(l) 
gauge symmetries, the electric charge of the resulting Dirac fermions is derived from a SU(2)vy 
quantum number (either +1/2 or —1/2 for a doublet) and the hypercharge Y as 



Q = Y + I 



3W, 



where I%w denotes the third component of "weak isospin" . Hence we can recognize the "up-" 
and "down-type" (anti)quarks of charge ±2/3 and +1/3 respectively, and (anti)leptons and 
neutrinos of charge +1 and 0. 



2.1 Higgs mechanism 

The reason why two weakly coupled gauge groups are needed rather than just the one of elec- 
tromagnetism is the existence of the weak interactions. The Fermi theory of weak interactions 
provides an adequate description of low-energy decays. But it is non-renormalizable and has 
many free parameters, which means that its predictive power is strongly limited; moreover, the 
theory is actually inconsistent at high energy, since the lowest-order contributions to scatter- 
ing cross-sections exceed theoretical bounds from unitarity. The introduction of heavy charged 
vector bosons W accounts for the weakness and short range of the interaction, and improves 
the theoretical behaviour of fermion scattering to some extent. But theories with massive vec- 
tors also in general have problems with unitarity and nonrenormalizability, and fail to describe 
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"neutral current" interactions in which (for instance) a charged lepton and a neutrino scatter 
without exchanging any charge. 

The solution to these problems was the construction of a renormalizable gauge theory of 
electroweak interactions based on the group SU(2)xU(l), which implies four gauge bosons, 
written as Wjf, W 3 and B^. The mass of weak vector bosons arises from the Higgs mechanism, 
i.e. spontaneous breaking of gauge symmetry by a scalar field. Since the Higgs scalar with 
quantum numbers (1,2) 1/2 




is not invariant under electroweak gauge transformation, if the scalar obtains a nonzero vacuum 
expectation value (v.e.v.) the gauge symmetry will be hidden. Effectively, the vacuum becomes 
a medium through which the gauge bosons cannot propagate as massless particles (as also in 
the theory of superconductivity). The Higgs' hypercharge is such that one component has zero 
electrical charge: this component may obtain a nonzero value, ((f) ) — v/y/2, without affecting 
the propagation of the photon, which remains masslessj^] 

Starting from the non-Abelian gauge theory with a scalar v.e.v. one can rewrite the action 
of gauge bosons in unitary gauge where the three components of 4> that leave |(</>)| invariant, 
the so-called "Goldstone bosons" , are absorbed by the W and one combination of the W 3 
and B. These vectors obtain masses proportional to (<j>) and to the gauge couplings, while the 
orthogonal combination of W 3 and B^ remains massless: 

Z p = cos 9 w Wl + sin BwB^ 
= - sin 9wW 3 + cos OwB^ 

where 

tan#vK = — , e = gsinOw, Mw = 7;gv, Mz — — -77—, Gf — (VZv 2 ) -1 . 
g 2 cos t>w 

Here g and g' are the coupling strengths of the underlying SU(2) and U(l) respectively, while 
e is the (magnitude of the) electronic charge and Gf the Fermi coupling. The components of <f> 
are accounted for as follows: one excitation about the vacuum value is inherently massive and 
becomes the Higgs boson, a neutral massive scalar. Excitations along the other three directions 
in field space correspond to gauge degrees of freedom along which the potential V(<f>) is constant, 
and would be massless were it not for the presence of the gauge bosons. In unitary gauge these 
excitations become the longitudinal modes of the gauge bosons, which thereby obtain three 
degrees of freedom (hclicity ±1,0) rather than the two (helicity ±1) of a massless vector. 



2.2 Gauge couplings and flavour sector 

Experimentally the parameters of the electroweak theory are found to be a ~ 1/137, sin 2 0w ~ 
0.23, v ~ 247 GeV, from which the underlying gauge couplings and scalar v.e.v. may be de- 
duced. This is a significant simplification over the old Fermi theory, which allowed a large 
number of undetermined couplings. But we have not yet considered how the Standard Model 
fermions obtain mass. As written in Eq. (JT|) they cannot become massive because two different 
representations cannot pair up in a gauge invariant manner, except for the singlet neutrino 
which can be given a Majorana mass. However, the Higgs solves this problem since its quantum 
numbers allow gauge invariant "Yukawa couplings" to be written down 

VVukawa ~ -htj(<f> ■ h)e C j + h.c. H 

2 In fact, whatever nonzero value of ((f) we take, one can find one U(l) subgroup that leaves it 
invariant. A redefinition of fields then brings the situation back to the standard one with electric 
charge defined as usual. 
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where i, j = 1 ... 3 run over the fermion family number, and (a • b) is an SU(2) invariant 
product. I have only written the lepton term explicitly. With a nonzero Higgs v.e.v. we obtain 
effective Dirac mass terms 



where m e = (<f>)h. e , the couplings h and masses to being 3-by-3 matrices in family space. 

Considering the up- and down-type quarks and charged leptons, one could a priori have 27 
complex couplings, greatly increasing the number of parameters^Most of these parameters are 
not physical: by field redefinitions under a global U(3) 5 flavour symmetry which acts on each 
Weyl fermion representation separately, one can bring the charged lepton mass matrix and one 
of the quark mass matrices (say m") to diagonal real form. The remaining quark mass matrix 
will still in general be complex and non-diagonal. 

Then a further field redefinition which treats the up- and down-type parts of the SU(2) quark 
doublet differently, and thus is not a symmetry of the weak gauge interaction, brings all matrices 
into real diagonal form. This is done at the cost of introducing a 3-by-3 unitary mixing matrix 
into the charged current (W^) interactions between quarks, the Cabbibo-Kobayashi-Maskawa 
(CKM) matrix. This has 4 parameters: the 9 parameters of a general unitary matrix, reduced 
by 5 through redefining the Dirac fermions by complex phases. 

Then the flavour sector contributes the following parameters: three charged lepton masses 
m e , rrifi, m T ; down-quark masses m^, m s , m&; up-quark masses m u , n\., m t ; and the 3 quark 
mixing angles #1,2,3 and one complex phase Sqkm of the CKM matrixrl 

At this point we have defined the Standard Model so far as it can be treated by perturbative 
quantum field theory with weak coupling strength and at lowest order. Higher order corrections 
do not change the basic picture of electroweak theory. However, the strong interactions cannot 
be described in this way. In particular, at low energies they are dominated by nonperturbative 
effects that introduce new mass scales, in addition to the Higgs v.e.v. and SM fermion masses. 
Also, in order to consider more fundamental theories that could lead to functional relations 
between the SM parameters, we will have to calculate with exponentially large ratios of mass 
or energy scales, for which lowest-order calculations are completely inadequate. Both these 
problems will be addressed in the next section. 



3 Energy scale dependence of couplings and unification 
3.1 Loop diagrams and RG evolution 

In quantum field theory (QFT) calculations of scattering processes and bound state energies, 
perturbative contributions are ordered by powers of fi. Higher order contributions include higher 
powers of coupling constants which we denote as 3, and some integrals over four- momenta of 
internal lines in Feynman diagrams, representing virtual fields which do not appear in the initial 
or final state. Such integrals can be schematically written as 



where s, t, u are the Mandelstam variables (certain finite combinations of initial- and final-state 
four-momenta), to is a particle mass and / is some known algebraic function. Depending on the 
theory, the integrals may be finite, or may diverge as a logarithm or power of the upper limit 
of integration. Thus in general the integral must be regularized, or cut off. This occurs either 
by a mathematical procedure in which an arbitrary energy scale A is introduced by hand; or by 
a physical cut-off where the structure of the theory changes at a given, large energy scale M uv 
and it is not appropriate to consider point-like virtual particles above this energy. The virtual 

3 The neutrino sector introduces still more, in general both Dirac and Majorana masses. 

4 Nonzero neutrino masses typically contribute at least three more masses, three mixing angles and 
one complex phase, however this count is model-dependent. 
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6 



Will be inserted by the editor 



particles then have some tightly-bound substructure which becomes evident at a length scale 
1/M UV ; below this length one should start to consider the individual constituents as propagating 
fields instead. 

The first possibility, regularization by hand, can only be useful if the physical quantities 
that emerge from summing contributions at different powers of fi can be made independent 
of the arbitrary regularization scale to sufficiently good precision. Remarkably, in a class of 
QFTs including the Standard Model, this is the case: such theories are called renormalizable 
2 . A renormalized coupling g^ may be defined, which has a functional dependence on A 
such that measurable quantities calculated to any given order in gn no longer depend on the 
regularization. 

However, the result may include logarithms of E/m, where E is a typical energy scale of 
the process (controlling the magnitude of s, t, u). In order for perturbation theory to work, the 
size of terms should decrease sufficiently quickly as one goes to higher orders of the expansion 
in h and g 2 : hence the presence of large logarithms ln(E/m) may be problematic. 

The appropriate solution is renormalization group evolution of coupling constants. Essen- 
tially we define a renormalized coupling g^ depending on an arbitrary energy scale fi which is 
allowed to varyj^] In perturbation theory we may then calculate the value of at a different 
scale to any given order as long as ln(/2'//i) is small. Now taking ln(/i'//2) — > a differential 
equation is obtained 

where the "beta-function" on the RHS may be expanded in powers of g^ and will depend on 
the gauge group and matter content of the theory. For some process of interest occurring at an 
energy or momentum scale ft, the size of loop corrections is minimized, and thus perturbation 
theory is most accurate, by expanding in powers of g^: this value is obtained from the value g' 
at another scale (characterizing some other physical process) via the RG equation (pi . 

In QCD (SU(3) Yang-Mills theory with fundamental Dirac fermion representations) we have 
at first order 

«*)-^(t-¥) < 3 > 

where n/ is the number of quark flavours contributing, i.e. quarks with mass m q < firj The 
negative sign of [3(g) has two notable consequences. First, at large momentum transfer the 
coupling decreases monotonically and loop corrections become less and less significant. QCD is 
asymptotically free and perturbative at high energies: experimentally, the "running fine struc- 
ture constant" a s (Mz) = 9s(Mz) 2 jAn ~ 0.118. Second, on integrating Eq. [2] a constant A c 
with dimensions of energy enters: the solution of the one-loop RG equation is 

,_v _ 12tt 

a 'W-(33-2 n/ )ln(Ai a M2) ( j 

and we see that a s diverges as the energy scale approaches A c from above. The coupling be- 
comes strong, higher-order terms should be considered and non-perturbative effects eventually 
dominate, such as confinement and the formation of a chiral quark condensate, which cannot 
be described within this framework: the degrees of freedom at low energies are no longer gluons 
and quarks. 

Nevertheless the behaviour of QCD (up to the effects of quark masses) is determined com- 
pletely by specifying the value of A c , and all energy scales associated with strong coupling 
effects such as the formation of hadrons are then proportional to A c . Thus we obtain some 
handle on the functional dependence of strong-interaction physics on the value of a s defined at 
some high energy scale. The value of A c (defined in some particular renormalization scheme) is 
then found by comparing perturbative calculations to scattering data and is of order 250 MeV. 



5 The bar distinguishes it from the proton-electron mass ratio fi. 

6 Particles heavier than the RG scale contribute negligibly as their propagators are suppressed by 
powers of mass. 
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Fig. 1. Experimental running of a s = g^l^n with energy (PDG review, 2006) 
The SU(2)vi/ and U(l)y gauge couplings of the Standard Model also show scale dependence: 



where n g counts generations of fermions. We may derive the running electromagnetic coupling 
ttem(A) and weak mixing angle 0\y(p) (the quoted value 0.23 applies for fi = Mz). However, the 
fine structure constant a is defined in the limit of zero momentum transfer (indeed at momenta 
lower than m e c the beta-function vanishes) therefore possible cosmological variations of a are 
independent of the question of RG scale dependence. 

3.2 Gauge unification and SUSY 

It was noticed some time ago that the weakest Standard Model gauge coupling becomes stronger 
with increasing energy scale, while the strong coupling becomes weaker; one may naturally ask 
whether at some high energy Mx the three forces could have a single common coupling strength 
and a common "unified" identity as aspects of a single more fundamental interaction. 

The Standard Model already contains examples of broken or hidden symmetries. The 
and Z° bosons result from a partially broken gauge symmetry in the Higgs phase, where the 
nontrivial properties of the vacuum state hide the gauge invariance. The strong interaction, 
meanwhile, exhibits confinement, meaning that particles charged under the QCD gauge sym- 
metry can never be observed: the energy cost of isolating a colour-charged particle would be 
larger than that of creating new particles to neutralize its charge. These two possibilities may 
also occur in new physics: there may be a larger gauge symmetry that is spontaneously broken, 
or the known particles may actually be composite bound states of some theory formulated at 
higher scale. I will focus on the former possibility simply because it is easier to calculate the con- 
sequences, as the coupling strength of such a theory turns out to be small and nonperturbative 
effects do not directly contribute. 

The first step towards unification concerned the matter representations of the SM. It was 
proposed [5] that the charged lepton SV(2) W doublets and singlets should be unified with the 
quark doublets and singlets in fundamental representations of SU(4). Moreover, an SU(2)# 
symmetry, also broken at high scale, unites the up-(neutrino-)type and down- (lepton-) type 
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Fig. 2. RG dependence of gauge couplings in SM and SUSY unification (M. Peskin) 

weak-singlet fermions. At leading order such Pati-Salam unification enforces relations between 
the masses of elementary fermions. The most convincing example is the heavy fermions of the 
third generation pQ, for which one obtains at the unification scale 



m T tan f3 = mi, tan [3 = m t — m Vr , /i = Mx 



up to small corrections, where m VT is the Dirac mass of the tau neutrino and tan/3 is the ratio 
of the v.e.v.'s of two electroweak Higgs doublets required by the theory. The observed difference 
between r and b quark masses is then accounted for by different RG running down from the 
unification scale, where the interactions of the quarks and leptons split off from one another. 

Most attention has focused on models that allow the fundamental forces to be unified 
(GUTs). The simplest mechanism for gauge unification is analogous to the symmetry-breaking 
of the Standard Model, but with a larger gauge group, in which some massive scalars receive 
nontrivial v.e.v.'s and their Goldstone modes are "eaten" by gauge bosons, leaving only the SM 
gauge bosons massless. Common symmetry breaking schemes start with a unified group SU(5) 
orSO(10). 

There are two notable consequences for the SM parameters. First, there may be algebraic 
relations between the Yukawa couplings and hence the masses and mixings of the SM fermions, 
as in Pati-Salam unification. Such relations depend strongly on the sector responsible for gauge 
symmetry breaking. Second and more importantly for us, the effective "fine structure constants" 
of the SM groups at the unified scale Mx obey a fixed relation 

g 2 s {M x ) = 9 2 {M X ) = gl{M x ) = 9x{M x ), g\ = ^-g' 2 . 

Note here that while the coupling strength of non-Abelian groups arise directly from the self- 
interaction of vector bosons, the coupling strength of U(l) is known only up to a constant 
factor which also scales the charges of matter fields. One needs to identify the charges of unified 
matter representations under the U(l) subgroup with the known SM charges to obtain a correct 
"hypercharge normalization" , resulting in the factor 5 /3 which occurs for many (though not all) 
unified models. The experimentally known coupling strengths at p, = Mz can be extrapolated 
to higher energies to test the unification relation. 



The result (Fig. 3.2 left-hand graph) is not very accurate for the simplest unification 
models)^] Non-minimal models including extra, heavy matter multiplets or multiple stages of 
symmetry-breaking can repair the discrepancy, but lose predictivity for our purpose unless 
additional assumptions are made. 

However, another way of adding extra particles has the effect both of fixing grand unifica- 
tion unambiguously and improving other properties of unified theories. This is supersymmetry 
(SUSY), a symmetry that relates each fermion in a theory to a boson with the same charge and 
mass (and vice versa). Thus the total count of particles is doubled]^] The particle physics mo- 
tivation for introducing SUSY is the technical naturalness of the mass parameter in the Higgs 



7 Such models are ruled out also by predicting too rapid proton decay. 

8 It turns out not to be feasible for any SM fermion to be the superpartner of a SM boson. 
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potential V((f>). This was omitted from previous discussion since the parameters of low-energy 
physics can be found without knowing the exact shape of V(^). 

The scalar mass is the only quantity in the SM whose value receives quantum corrections 
that diverge quadratically with the cutoff. Considering a more fundamental theory defined at 
a scale Mx, successive contributions to the scalar mass-squared in the low-energy theory are 
of order M x times some factors of coupling strengths. But the Higgs mechanism cannot work 
as observed unless the mass parameter is of order (100 GeV) 2 . Thus the scalar mass in the 
fundamental theory at Mx must be fine-tuned to extremely high accuracy: moreover, the fine- 
tuning must account for the exact quantum corrections, which cannot be found in perturbation 
theory with sufficient accuracy. However, precisely these quantum contributions to scalar mass 
vanish within supersymmetry, due to cancellation between fermions and bosons. 

Since the particle content is altered by SUSY, so are the RG equations, and it turns out that 
the f it of the measured SM couplings with a single grand unified coupling is much improved 



(Fig. 3.2 right-hand graph) with a unification scale Mx — 2 x 10 16 GeV and a unified coupling 
ax = gx(Mx) 2 /4"7r ~ 1/24QNOW we do not observe superpartners with the same mass and 
charge as the SM particles: however, it is possible to break supersymmetry but still preserve 
both the technical naturalness of the Higgs sector and grand unification. Soft SUSY-breaking 
introduces masses of order 100 GeV ~ ((f)) for the superpartners. Intuitively, this change should 
not affect properties that depend on the very large ratio Mx / (<f>) ■ We will consider the possible 
relation between superpartner masses and the Higgs v.c.v. in the next section. 

Now we may invert Eq. Q with rif = 3 and differentiate to obtain the dependence of A c 
(and thus the mass scale of hadrons) on the unified coupling: 

^ ln T7^ = ^ A\na x + —A\n c — | - + -A\n — ^ (5 

M x 9a x 27 M% 9 M% w 

where the last two threshold terms on the RHS arise from the change in the beta-function when 
the RG scale p, passes through the mass of a charged particle. Here rriq, rrig are the masses of 
the quark and gluon superpartners. If we similarly calculate the effect of RG running from Mx 
down to zero momentum on the electromagnetic coupling we find 

A\na = —Alnax+aY 9lIiA]n^- (6) 
3a x ^ 2tt M x W 

where the threshold term runs over all electrically charged fields with masses TOj. Note that in 
both cases the threshold terms are higher order in a or ax- Hence formally to leading order, 
the proton mass variation is expected to be [5] 

m 7T 

Aln— v - ~ — Z\lna~36zHna. (7) 
M x 12a v ' 

For instance, if there is a nonzero variation of 5 x 10~ 6 in a at some epoch, the corresponding 
variation of m p /Mx is about 1.8 x 10 -4 , which lies well outside observational bounds on /i = 
m p /m e variation from spectroscopy [B] unless the variation of the electron mass Aln(m e /mx) = 
(35 ± 5)Z\lna. We will however see that accounting for possible variations of SM fermion and 
superpartner masses can radically change the relation between variations of fi and a. 



3.3 Hierarchy problem and varying mass scales 

So far the origins of the Higgs v.e.v. ((f)) /Mx — 10~ 14 and small Yukawa couplings (e.g. 
rrid/rnb — 1/700) have not been discussed. The first of these small ratios is known as the 
hierarchy problem. As already discussed, SUSY stabilizes the hierarchy by ensuring the Higgs 
potential receives only logarithmic corrections proportional to the superpartner masses. An 



The unification was probed at LEP years after the first SUSY-GUT models were constructed. 
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alternative route is for the Higgs to be a composite state of a theory that becomes strongly 
coupled around the electroweak scale, A^c ~ 1 TeV, analogously to the strong coupling regime 
of QCD. In such technicolor models the theory above Atc no longer includes the Higgs and, as 
a weakly-coupled gauge theory, has only a logarithmic dependence on physics at higher energy 
scales. 

Then the ratio (</>) /Mx will depend exponentially on the technicolor gauge coupling at high 
energies: 

■M- = k a n h e-^ m ' bhah{M ^ (8) 
Mx 



where bh is the first term in the expansion of the beta-function (Eq. (|2])) with a factor — g 3 /167r 2 
removed: bh > means asymptotically free, cf. Eq. ([3]). Then if k and n are order 1, we may 
estimate m/bhOt^Mx) — 32 independently of the exact model used |7l8j given that Mx is 
order 10 16 GeV. Hence Aln{<p)/M x ^ 32Alna h (M x ). 

Within SUSY, essentially similar models are considered for a "hidden sector" where SUSY is 
broken at a scale exponentially below M x , in order to produce superpartner masses m of order 
the electroweak scale. Since quantum contributions to the Higgs potential are comparable to 
the superpartner masses, fh should not be much larger than (</>) except in special cases. In many 
unified models electroweak symmetry-breaking is caused exactly by such " radiative corrections" 
and one may thus set Aln((f>) ~ /Alum. The relation of Eq. ^ applies here also where ah now 
applies to the SUSY-breaking hidden sector. If the hidden sector or technicolor group is then 
unified at the GUT scale we obtain a relation between the variations of the Higgs v.e.v. and 
superpartner masses, and that of the unified coupling. 

However, the scale of electroweak symmetry-breaking may be strongly dependent on details 
of the model, including the variation of the top Yukawa coupling [5] which is a priori unknown. 
Hence we will keep the variation of ((f)) and fh as adjustable parameters: 

_ Alnm/Mx) _ A\n(m/M x ) . . 

Pv = -r-, , Ps = -r~i , (9) 

Zi in ax A in ax 

bearing in mind that they are likely to be of order 30 in the simplest unification scenarios. We 
can insert these variations in the Higgs v.e.v. and superpartner masses into the variations of 
A c and a, Eqs. (5|6l, to obtain [8] 



Ahifi 0.54/qx - 0-6/3, + 0.35fe 
~ A In a ~ 0.022/a x + 0.002/3,, + 0.011/3 S 

where for non-supersymmetric theories we set /3s = 0. Here we have also taken account of the 
small contributions of light quark (it, d and s) masses to the proton mass, which mean that m p 
is not exactly proportional to A c . 

Now for f3 v = Ps = 0, we obtain R ~ 25, independent of ax- This is a smaller variation of 
fj, than Eq. ^ due to the more accurate treatment of the proton mass, but still observationally 
problematic. 

For P v ~ 32, Ps = 0, and the choice ax — 1/24, we obtain R ~ —11 with an uncertainty of 
a few. With ax — 1/42, which may correspond more closely to non-SUSY unification, we have 
R ~ 3.5 with an uncertainty of about 5. 

In the case of varying superpartner masses, P v = Ps — 32 and ax — 1/24, we find R ~ 5±4. 
Hence the relation between varying \i and a is strongly sensitive to possible variations in the 
hierarchy of the Higgs v.e.v. or superpartner masses relative to Mx- Only a relatively small 
range of values of P v and Ps is compatible with a nonzero variation of \x being the same size as 
a variation of a. 



3.4 String theory 



The relation between particle theory and string theory is complicated and cannot be summa- 
rized briefly. I will only make a few general remarks on the status of coupling strengths and 
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mass scales in string models. Contrary to the field theory of point particles, string theory has 
one dimensionful fundamental parameter a', the inverse string tension (i.e. energy per unit 
length) of dimension [-E] -2 . One can equivalently define the "string scale" Ms ~ 1/ya' '. 

At energies much below Ms, only the lowest energy string excitations contribute to the 
dynamics: these are expected to behave like point particles to good approximation, and one 
can construct an effective field theory which describes their interactions in a flat space-time 
background, by matching field theory scattering amplitudes to string theory amplitudes. This 
is the starting point for much string model-building. One notable feature at this level is that 
the string coupling constant — the amplitude for one string to split into two — is given by the 
value of an effective scalar field, the dilaton. Thus we would need to know the dilaton potential, 
or at least its minimum, to make any statement about couplings. 

Supersymmetric string theory is, furthermore, necessarily defined in 10-dimcnsional space- 
time. Hence it is necessary to make six spatial dimensions invisible to low-energy excitations. 
This is achieved by compactification: a six-dimensional space is found which, in conjunction 
with 4-d spacetime, satisfies the equations of motion, and only modes with zero momentum 
along the six extra dimensions have sufficiently small mass to be experimentally accessible 

The details of four-dimensional physics then depend strongly on the six-dimensional space 
(and any fields that have nonzero components along it). In many cases there are continuous 
families of such spaces, parameterised by degrees of freedom that correspond to more scalar 
fields in 4-d spacetime, the moduli. Then we are interested in the minima of the potential as a 
function of the dilaton and the moduli; or if there is no minimum (as for an exponential potential 
[TU]L the speed of possible cosmological evolution of the fields. Recent developments in string 
model building have shown that most such scalars can become massive, i.e. the potential has 
minima in these directions. But the stabilization of all moduli and their masses and couplings 
to SM fields is an open question. We will see in the next section that the couplings of any light, 
cosmologically evolving scalar to the SM must be very weak (even weaker than gravity); this 
seems unlikely to be satisfied in most string models. 

In principle in a given string model, once the potential of the moduli is known, algebraic 
relations between the SM gauge couplings and Yukawa couplings should follow. In practice no 
model combines (semi-)realistic gauge symmetry matter and content with stabilized moduli to 
the extent of making predictions. Most "semi- realistic" models resemble grand unified theories 
in that there are simple relations between gauge couplings at the string scale, which however 
generally depend on the values of moduli. The details of how the SM gauge groups emerge and 
how the Yukawa couplings are determined are, however, different in string models. The main 
advance of strings over GUTs is that string models explicitly include gravity: the ratio Mg/M-p 
(the strength of gravity) is predicted as a function of the dilaton and moduli, leading to an 
extra condition on the values they may take consistently with observation. 



4 Deviations from GR and cosmological variations 

In particle theory we are required, due to the many precise tests of Lorentz invariance in 
microscopic physics [TT] to start from an action that is Lorentz covariant. Each of the fields 
representing a particle is a representation of the Lorentz group and the fields may only be 
combined in certain ways. This ensures that physical results will be unaltered under rotations, 
translations and boosts, after having controlled for environmental factors that do not respect 
this symmetry. Conservation laws such as energy and momentum are thus also built in. 

A situation where fundamental parameters are not unchanged under translation in time 
or space is then at first sight problematic: what becomes of the symmetry and conservation 
principles that give physics much of its simplicity and predictivity? For instance, merely intro- 
ducing (in QED) a fine structure constant a(x) that depends on space-time position will break 
energy and momentum conservation and gauge invariance, and removes any predictivity from 
the theory, since the function a(x) is a priori completely undetermined. 

10 This means that the size of compact extra dimensions should be at least as small as TeV -1 , except 
in a subset of "brane world" models where the SM cannot propagate along some direction. 
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Instead, we replace the unknown function by a dynamical field, a new degree of freedom 
whose behaviour is determined by its own action and couplings to the rest of the theory. We 
have already seen examples of this in the Higgs mechanism, where the SM fermion masses are 
functions of the Higgs v.e.v., and in string theory, where every coupling strength is due to 
some background field taking a particular value. Considering the action as a whole, Lorentz 
covariance is still valid: the variation in measurable quantities arises because the solution for 
the extra degree of freedom has a nontrivial space-time dependence. I will refer to such a field as 
a cosmon and denote it generically as ip. (It is assumed to be a scalar since other types of field 
generally cause gross violations of Lorentz symmetry if given a vacuum value.) An apparent 
violation of energy conservation in the SM, for example, would actually be a transfer of energy 
to or from ip, and should be correlated with other (in principle) measurable effects. 

Our task is then to relate the action of ip to observable effects. Although the action contains 
a priori unknown functions of tp, it may have predictive power in that the same configuration 
of the scalar should always produce the same observable effect: we assume that observable 
quantities are smooth and well-behaved functions of ip. 



4.1 Equivalence principles and their violation 

Considering variations with respect to gravitational as well as microscopic physics, we have 
to allow for the observational success of GR |12ll3j on scales from the laboratory to beyond 
the Solar System, which strongly restricts the theories we may write down. Note again that 
the appearance of space-time dependent fundamental parameters violates most of the basic 
principles of the theory. Consider first the Einstein equivalence principle (EEP): it consists of 
three parts. 

1. Weak equivalence principle (WEP): The trajectory of a freely falling test body is indepen- 
dent of its structure and composition. 

2. Local Lorentz invariance (LLI): The outcomes of non-gravitational experiments are inde- 
pendent of the velocity of their (freely-falling) reference frame. 

3. Local position invariance (LPI): Outcomes of non-gravitational experiments are independent 
of their position in space-time. 

The mass-energy of a test body depends on the fundamental parameters, and it is found 
(Nordtvedt) that in any situation where these are varying over spacetime the body will un- 
dergo an acceleration in addition to the gravitational acceleration due to GR. In general the 
masses of different test bodies will have different dependences, thus WEP is violated. The 
gradient of any given parameter defines a direction in spacetime, thus the boost and rotation 
invariance implied in LLI is broken, as (obviously) is LPI. 

The strong equivalence principle (SEP) is analogous to EEP, except that the universality 
of free fall holds also bodies with non-negligible gravitational self-energy, and the invariance 
principles apply also to local gravitational experiments. Violation of SEP occurs specifically if 
Gn varies (in units where particle masses are constant), for example in Brans-Dicke theory. Then 
the acceleration of a falling body will have a contribution from the gradient of its gravitational 
self-energy, which will depend on its mass and shape. The results of gravitational experiments 
(but not necessarily others) are also clearly dependent on space-time position and velocity in 
this case. 

The anomalous accelerations that signal non-GR effects may be simply expressed in terms 
of the gradients of the mass-energy of freely falling particles: 

- M dln(M/Mp)^ , . 

<Sa = -Vln — = A '- P; VlnG fc 11 

M P alnCxfc 

where M is the total mass-energy including binding energy and self-energy, and Gk are the 
values of fundamental parameters on which M may depend. In what follows we will neglect the 
gravitational self-energy, since this component is negligibly small for bodies whose accelerations 
have been most accurately measured. 
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Already at this level one can compare two methods of obtaining limits on variations of a 
within the Solar System. First via the dependence of atomic clock frequencies [TJ] as the Earth 
moves through the gravitational potential in its elliptical orbit; or second, in an Eotvos type 
experiment, via the differential acceleration 77 = \Sai — 80,2] /\g\ of test bodies with different 
fractions of electromagnetic binding energy B em /M = — d In M/dot. We have (up to a sign) 

ain(Mi/M 2 )|Vlna| 

"=^In^ T~' { ) 

Note that the last factor on the RHS is the ratio of the gradient in a to that in the Newtonian 
potential. Now despite recent improvements in atomic clock techniques, the Eotvos method 
is currently far more sensitive to spatial variations [15] . The limits on 7/ for various pairs of 
test bodies are at the 1CT 13 level [E] with A 12 (B cm /M) of order lCT 3 -lCr 4 . For the elliptical 
orbit of the Earth, the variation in the gravitational potential is a few times 10~ 10 , hence the 
seasonal variation in a consistent with bounds on rj is at most of order 10 -19 , barring some 



unusual cancellation in Eq. (11) |17j . Given that the clock frequencies vary as u>i oc a™ with 
rii of order 1, frequency stability at the 10~ 19 level would be required for the clock method to 
be competitive. Alternatively clocks could be sent into deep space to take advantage of a much 
greater variation in the 

The absolute size of anomalous accelerations Sa is also limited by tests of planetary motion 
and the form of the metric within the Solar System. This bound is relevant to Brans-Dicke theory 
where all test bodies experience the same acceleration due to the scalar interaction: the inverse 
scalar coupling strength uj in the theory is limited to be greater than 40, 000 However, in 
theories where particle physics parameters such as a and [i vary, differential accelerations give 
a much stronger bound. 



4.2 WEP violation via scalar couplings 

To relate local bounds on rj to more fundamental theories and to cosmological models we should 
write an action to describe the dynamics of the scalar (p = ip/Mp (here made dimensionless) 
and its couplings. By means of metric and field redefinitions this can be put into a form with 
the standard Einstein- Hilbert action for gravity and a canonical kinetic term for the scalar: the 
only remaining freedom is in the coupling functions of the scalar to SM matter and gauge fields. 
We may write a low energy effective action (strictly, action density) appropriate for dynamics 
at large distance such as gravitation: 

M 2 

£=^(R + dpipWip) - M*V{(p) - V vm + C cm + ■■■ (13) 

where I have omitted kinetic terms for matter fields, and nuclear interactions. Here V^m contains 
the interactions of the scalar with particles that constitute ordinary matter, while C em gives 
the dependence of the electromagnetic interaction: 

V vm = TO e0 c e (^)ee + m N0 c N (<p)(hn + pp) + ^6 N0 c v (p)(nn - pp) (14) 
-£ em = ^Z{0)F^F" v , (15) 

where to^o = (jn n + m p )/2 and 6 no — rn n — m p for present-day values of masses, thus the 
functions Ci and Z are normalized to unity for the present value tpo- Then 

41na = -lnZ(£), Z\ln/^ = ln V ' , 2mwo K - (16) 

For a small change in (p about the present value we expand = 1 + Ai({p — <po) + ■ ■ ■ to define 
the couplings A^jy.,, and Z((p) — (1 + A cm (<p — i^o))^ 1 - [""[Thus Zilna = X cm A(p and Z\ln/i = 

11 Under some special assumptions [20], the leading term in the expansion of A; and Z~ x about some 
particular value may be quadratic in (p, then the scalar forces will be much suppressed. 
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(A, — X v 8jsr /2mjv — X e )Aip. Comparison with Brans-Dicke theory, which may be reformulated 
as a scalar with universal coupling (3 + 2cj) -1 / 2 to matter, indicates that the magnitude of 
couplings is bounded at a few times 10~ 3 before considering differential acceleration. 
We then consider the dependence of atomic mass per baryon M / A on the scalar: 

-j = f P (™ P + m e ) + (1 - f p )m n + - Bst ™ a z (17) 

where the proton fraction f p = Z/(Z + A). Here E Dm is the nuclear electromagnetic self- 
energy, estimated as (a/ao)f em mjs[Z(Z — l)/^! 1 / 3 with / om ~ 7 x 10~ 4 , and -B s tron g is the 
nuclear binding energy due to strong interactions. Owing to the complex and poorly known 
dependence of nuclear forces on fundamental parameters, this last contribution has generally 
been neglected. To find the differential acceleration of two test bodies towards a given massive 
source, the source coupling dln(M/AMp) s /d(p, which gives the value of through solving a 

Poisson equation for <^(r), and the test body couplings dln(M/AMp)i^/d<p, must be evaluated. 
The source coupling is always dominated by din cat / 'dip ~ Ajv 5 thus we may write 



~2~ 



(—A w h X e fX 1 )Ai2.f„ + A cm / cm Z\i2 



(18) 



where A\% denotes the difference between test bodies: typically Anfp may be of order 0.03 
while A 12 (Z 2 /A 4 / 3 ) is of order lpl 

Within any given unified model, we may define Xx = din ax /dtp and calculate the low- 
energy couplings using previous results for nucleon and electron masses and a |8] : the result for 
rj is then proportional to X 2 X [3T]. Next we relate this bound to cosmological and present-day 
variations of a and p. 



4.3 Scalar evolution in cosmology and space-time variation 

Clearly by making Ax arbitrarily small any bound from differential accelerations may be evaded; 
but the cosmological variation also becomes negligibly small, unless the variation of (p is suf- 
ficiently large. Here we consider the cosmological effects of an evolving scalar. While the local 
gravitational effects may be estimated by solving for the scalar profile in a static spacetimc, 
large-scale cosmological effects are found by considering the time evolution of an approximately 
homogeneous and isotropic metric of Robertson- Walker form ds 2 = dt 2 — a 2 (t) dx 2 . The 
equation of motion of the scalar is then 

(p + 3Hp = -V'(<p) - Kff(^ Pmattor) (19) 

where H = a/a is the Hubble expansion rate. Here V^g(ip,p m ) ~ p m d\sxm/d(p arises from 
the coupling of the scalar to (nonrelativistic) matter, including possible dark matter, m being 
the mass of the relevant particle. If we assume small couplings, the bare potential V'(ip) will 
generally dominate. Note that the effective mass-squared V"(ip) should be extremely small 
(comparable with H 2 ) if the value of the scalar is to drift monotonically over cosmological time 
rather than oscillate about a minimum. The kinetic energy of the field is T — <p 2 /2, resulting 
in an effective equation of state w v = p v /p v = (T — V)/(T + V). 

The evolution of a scalar may be connected with the most conspicuous apparent deviation 
from GR on large scales: this is the recent accelerated expansion of the Universe, a > 0, required 
to explain a number of independent cosmological observations: see e.g. [55]. Acceleration re- 
quires a total equation of state w < —1/3 averaged over all mass-energy components, whereas 
matter (baryonic or dark) has w rn ~ and radiation has u> 7 = 1/3. Thus some additional 
form of energy density with negative pressure should be present and dominate over matter and 
radiation in recent cosmological evolution. This is the so-called "dark energy" . 

12 This treatment is valid provided that the scalar's Compton wavelength (inverse mass) is larger than 
the distance between source and test bodies. Otherwise we have a suppression factor e - ™* 5 . 
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Fig. 3. Observational bounds on matter density and dark energy equation of state, from |23] . A 
spatially flat Universe is assumed, for which J?de + &m — 1-0. 



The theoretically simplest candidate for dark energy is a cosmological constant A, which 
was already allowed by Einstein in the formulation of GR. The corresponding equation of state 
is u>a — — 1, which is compatible with all current observations. However, it is also possible that 
scalar fields have contributed significantly to the energy density either at present or at earlier 
epochs. We may either assume that the entire acceleration is due to ip, or use observations to 
put bounds on the rate of change <p at recent epochs. 

The rate of variation is related to the energy density fraction Q v = 



cp _ V 21 
H = HM P 



yj3n v (l + w v ) (20) 



where both quantities inside the square root are bounded above observationally. Recent mea- 
surements give an allowed region in the Q v -w v plane, Fig. [3] within which Eq. (20 1 attains a 



maximum value of approximately 0.7. Thus the recent scalar evolution cannot have been faster 
than <p max ~ 5 x 10~ n y -1 , using Hq ~ 7 x I0~ 11 y~ 1 . This bound was obtained assuming 
that <p accounts for the entire dark energy, but it holds much more widely because 1) the total 
magnitude of dark energy density J?<ie is the relevant quantity affecting cosmological evolution, 
and bounds fl v above; and 2) the scalar kinetic energy influences the total equation of state 
whether or not ip constitutes the dark energy. As an extreme example, a scalar mimicking cold 
dark matter with w v ~ and f2 v ~ 0.3 would have <p>/H ~ 1. 

We may now write dthi{a, fi} < {ci, C2}Xx ^max where C\ t 2 are numerical constants that 
depend on the choice of unified model. For the differential acceleration we have 

V^K^X% (21) 



where K is another model-dependent constant arising from Eq. (18 1, where we keep only the 
first term on the RHS which dominates for the models we consider. Thus for a nonzero variation 
of a or fi, setting Ai 2 f p — 0.036 we find 

K( ft In a \ 2 K ( d t ln^ \ 2 

V - c\ ^3.7x10-1° y~V ' V - c| \3.7x 10"i° y"V ' (22) 

For instance, in a scenario where the unified coupling varies but the Higgs v.e.v. and superpart- 
ner masses do not, V — /3s = in Eq. d9|), we have — {0.52, 13} and K ~ IX^XydN/mN] — 
0.18, thus K/ci 2 — {0.67,1.1 x 10 _3 }T3ence a variation a/a ~ 10^ 16 per year would lead 
to rj < 5 x 10 -14 . Eotvos experiments may be competitive with atomic clocks in bounding 
variations in such models |21l24j . 

Conversely, in a scenario where where (3 V ~ (3s — 34, we find — {1-0,4} and K ~ 1.5, 
thus K/c\ 2 = {1.5,9.5 x 10~ 2 }. The small variation of /i is due to a cancellation between 
variations of m p /Mx and m e /Mx, whereas the scalar couplings to the proton and electron 
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remain relatively large. Thus for a linear time variation fi/fj, ~ 2 x 10 y , which would 
yield Afi/n ~ 2 x 10~ 5 at high redshift, this scenario predicts rj > 10~ 12 , already beyond 
experimental bounds. 

Several of the assumptions we made in deriving these results may be violated: for example 
scalars may be significantly heavier than the Hubble scale, or may have couplings to matter 
or self-couplings that are not weak. The effective scalar mass may well depend on p m , which 
is the special feature of so-called "chameleon" models [25 whose gravitational effects are quite 
different from very light, weakly-coupled fields. In such cases the dominant variation of the 
scalar is likely to be with respect to spatial position or local density [26] rather than a slow 
cosmological drift. The framework of an effective action for the scalar and matter fields may also 
be used in these cases to predict the astrophysical or cosmological behaviour of such theories. 



5 Conclusion 

I have sketched the theoretical framework of the Standard Model and General Relativity nec- 
essary to understand the role of "fundamental constants" in theoretical physics, and discussed 
theoretical extensions which could account for the values of these parameters and their possible 
cosmological variation. There are many potential and actual observational probes of such new 
physics associated with "varying constants" . Determining their history brings together many 
branches: astrophysics, optics, gravitation, and atomic, nuclear and particle physics. At present 
we are still awaiting a clear signal of non-standard behaviour; nevertheless the search for varia- 
tion, and the associated violation of equivalence principles, involves some of the most sensitive 
and precise experiments in physics. 
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